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ABSTRACT 


Recent  observational  data  cast  serious  doubt  on  the  widely  held 
view  that  the  sun's  corona  is  heated  by  traveling  waves  (acoustic  or 
magnetohydrodynamic).  It  is  here  proposed  that  the  energy  responsible 
for  heating  the  corona  is  derived  from  the  free  energy  of  the  coronal 
magnetic  field  derived  from  motion  of  the  "feet"  of  magnetic  field  lines 
in  the  photosphere.  Stochastic  motion  of  the  feet  of  magnetic  field  lines 
leads,  on  the  average,  to  a  linear  increase  of  magnetic  free  energy  with 
time.  This  rate  of  energy  input  is  calculated  for  a  simple  model  of  a 
single  thin  flux  tube.  The  model  appears  to  agree  well  with  observational 
data  if  the  magnetic  flux  originates  in  small  regions  of  high  magnetic 
field  strength  as  proposed  by  Tarbell,  Title  and  Schoolman.  On  combining 
this  energy  input  with  estimates  of  energy  loss  by  radiation  and  of  energy 
redistribution  by  thermal  conduction,  we  obtain  scaling  laws  for  density 
and  temperature  in  terms  of  length  and  coronal  magnetic  field  strength. 


I.  INTRODUCTION 


Following  the  general  acceptance  of  the  concept  of  a  very  hot  corona 
in  1945  (Billings,  1966),  Biermann  (1948)  and  Schwarzschild  (1948)  inde¬ 
pendently  suggested  that  the  corona  is  heated  by  acoustic  waves.  However, 
recent  observational  data  appear  to  make  this  hypothesis  untenable. 

Athay  and  White  (1979),  in  their  analysis  of  UV  spectroscopic  data 

obtained  by  means  of  the  0S0-8  Spacecraft,  argue  that  the  acoustic  wave 
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flux  in  the  chromosphere  cannot  exceed  about  10  erg  cm  s  whereas  the 
radiation  losses  of  the  transition  region  and  corona  require  an  energy 
flux  of  10  '  erg  cm  s  .  Bruner  (1980),  from  an  independent  analysis 
of  0S0-8  UV  spectroscopic  data,  argues  that  most  of  the  wave  motion  in  the 
transition  region  is  in  the  form  of  standing  waves  (probably  evanescent 
waves)  rather  than  traveling  waves,  so  that  the  net  acoustic  flux  is  at 
least  three  orders  of  magnitude  lower  than  that  needed  to  heat  the  corona. 

The  presence  of  a  magnetic  field  complicates  the  discussion.  It  is 
possible  that  the  change  is  minor,  in  that  energy  still  propagetes  as 
waves,  but  these  are  magneto-acoustic  waves  and/or  Alfven  waves  (Osterbrock, 
1961;  Stein  and  Leibacher,  1974).  If  the  waves  are  magnetoacoustic,  we 
still  have  a  problem  in  accounting  for  the  required  energy  flux.  If  the 
waves  are  Alfven  waves,  then  there  is  a  problem  in  accounting  for  their 
dissipation  (Stein  and  Leibacher,  1974)  although  Uchida  and  Kaburaki  (1974) 
have  argued  that  large-amplitude  Alfven  waves  may  convert  to  magneto¬ 
acoustic  waves  at  coronal  heights  and  thereby  be  dissipated. 

On  the  other  hand,  the  close  association  between  magnetic  field 
strength  and  coronal  heating  suggests  that  the  magnetic  field  may  play 
an  active  role,  rather  than  a  passive  role,  in  the  energy  transport 
process.  For  instance,  it  is  notable  that  the  Rosner-Tucker-Vaiana  model 


of  hot  coronal  loops  (Rosner,  Tucker  and  Vaiana,  1978)  seems  to  fit  not 
only  large-scale  loop  structures  and  active  region  loops,  but  also  coronal 
loops  produced  by  solar  flares.  Although  this  model  deals  only  with  the 
relationship  between  energy  transport  by  heat  conduction  and  by  radiation, 
it  does  raise  the  question  as  to  whether  coronal  loops  are  heated  by  some 
mechanism  akin  to  that  responsible  for  solar  flares,  which  is  believed  to 
be  the  release  of  free  magnetic  energy  by  dissipative  plasma  processes 
(Sturrock,  1980).  It  is  quite  possible  that  the  explosive  nature  of  flares 
is  due  to  the  sudden  rearrangement  of  magnetic  field  by  an  MHD  instability: 
heating  (and  possibly  acceleration)  which  occur  during  a  flare  may  be  due 
to  dissipative  processes  which  may  also  occur  at  a  slower  rate  in  the 
"quiet"  solar  atmosphere. 

It  appears  that  there  are  two  requirements  for  the  propagation  of 
preflare  energy  from  the  photosphere  into  the  corona:  there  must  be  a 
magnetic  field,  and  the  footpoints  of  the  magnetic  field  must  be  moved  in 
such  a  way  that  the  configuration  is  raised  in  energy  from  the  current- 
free  "ground  state"  to  the  current-carrying  "excited  state"  which,  in  the 
early  stage  of  excitation,  may  with  good  approximation  be  taken  to  be  a 
force-free  magnetic-field  configuration.  These  requirements  can  be  met 
equally  well  in  the  "quiet"  solar  atmosphere:  the  photosphere  is  always 
permeated  by  magnetic  field,  and  the  photosphere  is  always  in  motion  due 
to  granulation,  supergranulation  and  other  motions.  The  aim  of  this 
article  is  to  estimate  the  resulting  flux  of  energy  from  the  photosphere 
into  the  corona  and  to  incorporate  this  assumption  into  the  Rosner-Tucker- 
Vaiana  model  so  as  to  obtain  estimates  of  the  density  and  temperature  of 
the  plasma  in  a  coronal  loop  in  terms  of  the  basic  parameters  of  the  loop. 

This  possibility  has  recently  been  addressed  also  by  Golub  et  al . 
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(1980).  They  begin  with  an  analysis  of  observational  data  which  indicates 
that  there  is  a  relationship  between  the  plasma  pressure  in  a  coronal  loop 
and  the  magnetic  field  strength  in  the  loop,  which  they  assert  to  be 
p  «  B^‘5.  They  estimate  the  rate  of  energy  supply  into  a  coronal  loop  in 
terms  of  the  mean  transverse  (toroidal)  magnetic  field  strength  and  the 
mean  transverse  (torsional)  velocity  at  the  photosphere,  but  without  con¬ 
sideration  of  the  cause-and-effect  relationship  of  these  two  quantities 
or  their  stochastic  nature.  Analysis  of  the  data  contained  in  their 
Fig.  2(b)  yields  a  pressure-magnetic-field-strength  relationship  of  the 
form 

p  «  Bb,  b  =  0.77  ±  0.23.  (1.1) 

This  is  a  better  fit  to  the  relationship  obtained  in  this  article  (b  =  6/7, 
see  Section  5),  than  the  value  b  =  72/7  found  by  Golub  et  a  1.  (19 80). 
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II.  TWISTED  FLUX  TUBE 


We  wish  to  consider  the  response  of  a  coronal  flux  tube  to  photos- 
pheric  motion.  Since  the  free  energy  available  for  heating  the  plasma  is 
due  to  currents,  we  consider  in  particular  a  twisting  motion  such  as 
would  result  from  a  rotation  of  either  or  both  ends  of  the  flux  tube  at 
the  photosphere.  For  simplicity,  we  consider  a  tube  which  is  thin  compared 
with  its  length.  We  also  simplify  the  calculation  by  considering  a  tube 
for  which  the  central  field  line  is  a  straight  line  and  the  tube  has 
cylindrical  symmetry.  As  indicated  in  Section  I,  we  assume  that  the 
plasma  density  is  sufficiently  low  that  the  field  is  approximately  in  a 
force- free  state. 

It  will  be  seen  that  the  model  we  are  considering  is  very  close  to 
that  of  Gold  and  Hoyle  (1960)  except  that  we  are  allowinq  for  slow  vari¬ 
ation  of  the  field  with  respect  to  the  z  coordinate  of  cylindrical  coordi¬ 
nates  (z,r,4>).  On  the  other  hand,  we  will  simplify  the  model  by  consider¬ 
ing  only  the  lowest  order  significant  terms  in  polynomial  expressions  of 
quantities  in  terms  of  r. 

With  these  restrictions,  we  find  that  the  solution  of  the  equations 


7.B  =  0 


Bx(Tx  B)  a  0 


may  be  expressed,  to  lowest  significant  order  in  r,  as 


Bz(z,r)  *  Bq(z) 


-[iBo  (I)  * 


Br(z,r)  «  -  j  Bg  (z)  r, 


B.  (z,r)  *  b  Bq  ( 2 )  r. 


We  see  that  each  field  line  is  rotated  about  the  axis  by  an  angle  Ax 
given  by 

AX  =  bL,  (2.4) 


when  L  (cm)  is  the  length  of  the  tube. 

If  the  radius  of  the  flux  tube  is  R(z),  the  magnetic  flux  passing 
through  the  tube  is  given  by 

*  =  "Bo  r2  -  g-  B;1  R4  -  J  b2B0R4,  (2.5) 


and  the  magnetic  pressure  at  the  surface  of  the  tube  is  given  by 


d  -  1  R2  +  1 

P  S7  8o  87 


1/bM2 
14  (V 


1  BX' 

2  0  0 


(2.6) 


Let  us  now  consider  a  flux  tube,  which  is  initially  untwisted  so 
that  b  =  0,  immersed  in  a  plasma  so  that  the  magnetic  pressure  at  the 
surface  of  flux  tube  is  balanced  by  an  equal  gas  pressure  outside  the  tube. 
Now  suppose  that  the  tube  is  twisted  so  that  b  f  0,  but  we  require  that  the 
flux  $  is  unchanged  and  that  the  magnetic  pressure  p  at  r  =  R  is  unchanged, 
since  it  is  still  balanced  by  the  same  external  gas  pressure.  Then  we 
must  expect  that  BQ(z)  and  R(z)  will  both  change.  It  turns  out,  however, 
that  to  lowest  order  (quadratic)  in  b,  R(z)  is  unchanged.  On  the  other 
hand,  BQ(z)  is  found  to  increase  quadratically  with  b  according  to 


The  "free  energy"  AW  of  the  flux  tube  is  the  magnetic  energy  of  the 
twisted  tube  less  the  magnetic  energy  of  the  unt'./isted  tube.  To  lowest 
order  in  b,  this  is  found  to  be 


(2.8) 
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which  shows  that  the  free  energy  is  distributed  uniformly  along  the  length 
of  the  twisted  flux  tube.  By  using  equations  (2.4)  and  (2.5),  we  may  re¬ 
express  (2.8)  as 


III.  STOCHASTIC  MOTION 


We  now  consider  the  twisting  of  a  flux  tube  due  to  stochastic  motions 
at  the  photosphere.  The  general  surface  motion  may  be  divided  into  two 
components,  according  to  the  equation 

v  =  -  vy  +  7x(An)  (3.1 ) 

where  n  is  the  unit  vector  normal  to  the  surface.  The  first  term  on  the 
right-hand  side  denotes  the  "curl -free"  component  and  the  second  term 
denotes  the  "vortical"  component. 

We  need  to  consider  the  transfer  of  energy  from  the  photosphere  into 
the  magnetic  field  due  to  a  random  (horizontal)  velocity  field  in  the 
photosphere.  A  general  treatment  will  be  published  at  a  later  date.  At 
this  time,  we  undertake  only  a  simplified  analysis. 

We  assume  that  the  velocity  field  is  statistically  homogenous  and  iso¬ 
topic,  and  that  half  of  the  energy  in  the  velocity  field  is  in  curl -free 
motion  and  half  in  vortical  motion.  Only  the  latter  component  contributes 
to  twisting  of  magnetic  field  lines.  We  therefore  consider  an  elementary 
contribution  to  the  twisting  by  considering  the  rotation  of  one  end  of  the 
flux  tube  where  it  meets  the  photosphere.  Another  contribution  to  the 
stored  energy  will  come  from  the  other  end  of  the  flux  tube. 

If  the  vorticity  at  the  center  of  the  flux  tube  where  it  meets  the 
photosphere  is  i>, 

vh  =  -r  (3.2) 

and  the  contribution  to  the  twisting  due  to  rotation  at  one  end  only  is 
given  by 

t 

A-<  -  /  *(t')dt\  (3.3) 

o 
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This  leads  to  the  rate  of  change  of'energy  given  by 

dW  ,  t2  /  (Ax)2 
dt  1  6tt2L  At  ' 

We  see  from  equation  (3.3)  that 

t+At  t+At 


(3.4) 


This  becomes 


<(  (Ax)2)  «  /df  /dt-  <*(t‘  )u(t"  )>  . 
t  t 

00 

<)ax2)^)  ~  At  j dx  R(t) 


(3.5) 


(3.6) 


if 


R(t)  =  ( t  )co(  t  +  x)^> 


(3.7) 


and  At  is  large  compared  with  the  range  of  x  over  which  R(x)  is  significant. 


Hence  if  the  "correlation  time"  xc  is  defined  by 


I- 


'dx  R(x)  =  R(0)x 


and  if  At  »  x  ,  equation  (3.6)  becomes 


(3.8) 


.  2  r(0)tc  -  2  <„*>V 


(3.9] 


If,  considering  only  one  end  of  the  flux  tube,  the  radius  of  the  flux 
tube  is  R*  and  the  field  strength  is  B*  at  the  photosphere,  then 

=  V  (3'10) 

2  2 

The  combination  (j/ >  R*  is  the  mean  square  vortical  velocity  at  the 


circumference  of  the  flux  tube,  which  is  twice  the  mean  square  vortical 
component  of  the  velocity  over  the  flux  tube.  Hence  it  is  the  same  as 


2 

the  mean  square  velocity  (  v  )  ,  including  both  vortical  and  curl -free 
components. 

If  we  now  consider  both  ends  of  the  flux  tube  and  assume  that  the 
field  strength  is  the  same  at  both  ends,  we  find  that  the  rate  of  increase 
of  energy  in  the  flux  tube  is  given  by 


dW 

dt 


rar<u  >v 


(3.11) 


Now  suppose  that  Rc  is  the  mean  radius  of  the  flux  tube  in  the  corona 
defined  by 

V  =  ttR2L  (3.12) 

c 


where  V  is  the  volume  of  the  flux  tube.  Suppose  also  that  the  mean  field 
strength  Bc  is  defined  by 

t  =  ttR2  Bc  .  (3.13) 

Then  if  £j  (erg  cm  s  )  is  the  mean  rate  of  energy  input  into  the  flux 
tube,  defined  by 

^  =  ejV,  (3.14) 

we  see  that 

£;=  K  Bc  L"2  (3.15) 

where 

K  =  B*<v2>T(;.  (3.16) 
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IV.  SCALING  LAWS  FOR  CORONAL  LOOPS 


We  now  consider  a  coronal  loop  and  attempt  to  obtain  scaling  laws 
for  the  coronal  density  and  temperature  by  investigating  the  energy 

_3 

balance.  We  denote  by  nQ  (cm  )  and  Tq  (K)  the  density  and  temperature  at 
the  top  of  the  loop.  For  simplicity,  we  consider  only  loops  which  are 
sufficiently  small  that  they  are  approximately  isobaric.  Then,  if  the 
mean  density  and  temperature  are  nc  and  T  , 


We  follow  Rosner,  Tucker  and  Vaiana  (1978)  in  approximating  the  radi- 

O  1 

ation  energy-loss  function  eR  (erg  cm  s  )  by 

•  r  "c  Tc/2  <4-2> 

wherein  F  =  10"^'®. 

If  energy  is  deposited  primarily  near  the  top  of  a  loop,  it  would  be 

carried  to  the  lower  regions  of  the  loop  by  heat  conduction.  If  we  denote 
-3  -1 

by  eH  (erg  cm  s  )  the  rate  at  which  energy  is  extracted  from  the  upper 
regions  of  the  loop  by  heat  conduction  for  transfer  to  the  lower  regions 
of  the  loop,  then 

£H  =  0  To/2  L'2*  (4.3) 

Noting  that  the  "temperature  scale  height"  is  of  order  1/2  L,  and  that 
the  thermal  conduction  coefficient  for  a  fully  ionized  plasma  (Spitzer, 
1962)  is  approximately  10~®  T^,  we  see  that  9  =  lO"^'4. 

If  we  now  assume  that 


T 

c 


=  6  T 


o’ 


(4.4) 
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we  may  obtain  a  relationship  between  n  and  T  by  assuming  that  e„  ~  zu. 

0  0  K  H 

This  relation  is  found  to  be 

To  =  a-'^4  r'/4  a-5/8  l'/2.  (4.5) 

Since  the  mean  rate  of  energy  input  must  be  balanced  by  the  mean 
rate  of  energy  output  (by  radiation),  we  obtain  what  should  be  a  more 
reliable  relationship  by  equating  the  expressions  in  equations  (3.15) 
and  (4.2).  This  leads  to  the  relation 

n2  T"1/2  =  r"1  KB  L"2  (4.6) 

c  c  c 

for  n„  and  T  . 
c  c 

Alternatively,  pursuing  the  assumption  that  eR  ~  eH,  we  may  set 
eH  ~  £j  and  so  obtain,  from  equation  (3.15)  and  (4.3),  the  expression 


V.  DISCUSSION 


Equation  (4.5)  has  the  same  form  as  a  relation  obtained  by  Rosner, 
Tucker  and  Vaiana  (1978).  The  more  detailed  model  analysis  of  Vesecky, 
Antiochos  and  Underwood  (1979)  leads  to  a  similar  relationship  with  the 
numerical  form 


To  = 


10 


-3 


,25  nl/2 
0 


1/2 


(5.1) 


We  find  that  this  agrees  with  equation  (4.5)  if  we  adopt  the  proposed 

-0  15 

values  of  0  and  r  and  adopt  5  51  10  '  ’  i.e.  Tc/Tq  *  0.7,  which  seems 

not  unreasonable. 

Equation  (4.7)  is  of  special  interest  in  that  it  leads  to  a  relation¬ 
ship  between  Tq  and  Bc>  independent  of  the  length  L.  Although  we  have 
many  magnetograph  maps  of  the  solar  disk,  unfortunately  we  do  not  know 
the  magnetic  field  strength  at  the  top  of  a  coronal  loop.  Golub  et  al. 
(1979)  have  recently  presented  data  for  the  coronal  temperatures  of  three 
typical  features:  x-ray  bright  points,  active  regions,  and  large-scale 
structures.  The  data  presented  by  Golub  et  al.  (1979)  suggests  the 
approximate  numerical  relationship 

T  *  105'9  1  0,1  B2/7.  (5.2) 

u 


In  order  for  equation  (4.7)  to  have  this  numerical  form,  we  require  that 
K  =  ID15’2  ~  We  see  from  equation  (3.16)  that  this  requires  the 

following  combination  of  parameters  at  the  photosphere: 


<v2> 


10 


16.3  *  0.4 


(5.3) 


Tarbell,  Title  and  Schoolman  (1979)  propose  that  most  of  the  magnetic 
field  in  the  solar  atmosphere  arises  from  small  knots  of  intense  field  of 
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strength  8*  =  1200  and  that  the  r.m.s.  velocity  field  in  the  photosphere 
5  -1 

is  about  10  cm  s  .  Hence  our  theory  leads  to  the  empirical  relation¬ 
ship  (5.2)  if  the  correlation  time  for  the  photospheric  velocity  field 
has  the  value  tc  «  103-3  ±  0,4  s,  i.e.  in  the  range  10  -  80  m.  The 
lower  limit  of  this  range  is  comparable  with  current  estimates  (8  m)  of 
the  mean  lifetime  of  granules  (Allen,  19 73). 

At  this  time,  observational  data  concerning  photospheric  motions 
and  the  photospheric  magnetic  field,  and  observational  data  concerning 
magnetic-field  strength  at  coronal  levels,  are  insufficiently  precise  to 
provide  a  definitive  check  of  the  expression  (4.7)  for  the  coronal  temper¬ 
ature  in  terms  of  the  coronal  magnetic  field  strength.  In  any  case, 
relations  derived  in  Sections  IV  and  V  depend  on  highly  idealized  assump¬ 
tions  concerning  the  structure  and  energy  balance  of  coronal  loops.  The 
present  estimates  have  been  made  simply  to  show  that  the  heating  mechanism 
proposed  in  this  article  is  not  obviously  ruled  out  by  observational  data. 

In  order  to  obtain  more  detailed  consequences  for  comparison  with  observa¬ 
tional  data,  we  intend  to  incorporate  the  heating  mechanism  discussed  in 
Section  III  into  the  model  analyzed  by  Vesecky,  Underwood  and  Antiochos 
(1979).  In  this  way  we  shall  avoid  certain  restrictive  assumptions,  such 
as  constant  pressure  and  the  assumption  that  eH  =  eR,  and  we  shall  be  able 
to  include  the  important  effect  of  magnetic-field  geometry. 

Nevertheless,  the  fact  that  the  models  of  Vesecky,  Antiochos  and 
Underwood  (1979)  lead  to  a  scaling  law  similar  to  that  derived,  on  a  simpler 
basis,  by  Rosner,  Tucker  and  Vaiana  (1978)  give  some  cause  for  optimism 
that  the  scaling  laws  derived  on  the  basis  of  the  restrictive  assumptions 
of  the  current  article  may  also  prove  to  have  wider  applicability  than 
might  initially  be  expected.  Ignoring  the  numerical  coefficients,  we 


see  from  equations  (4.6)  and  (4.7)  that  the  density  and  temperature  in 
a  coronal  loop  will  scale  with  the  length  of  the  loop  and  the  coronal 
magnetic-field  strength  as 


n  a  B4/7  L-1 
c  c 

(5.4) 

V  Bc/7' 

(5.5) 

(5.6) 

so  that 


It  is  notable  that  the  index  6/7  (0.86)  appearing  in  the  scaling  relation 
(5.6)  is  quite  consistent  with  the  value  0.77  t  0.23  derived  from  the  data 
of  Golub  et  al.  (1980),  as  noted  in  equation  (1.1),  and  is  in  fact  a 
better  fit  than  the  value  12/7  which  they  derived  from  their  theoretical 
analysis. 

It  will  be  interesting  to  see  whether  new  data,  such  as  may  be  obtained 
from  the  Solar  Maximum  Mission,  will  substantiate  the  above  scaling  laws. 
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